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ABSTRACT 

We explore the possibility that extended disks, such as that recently discovered in M31, are the 
result of a single dwarf (10 9 -10 10 A/ Q ) satellite merger. We conduct N-body simulations of dwarf 
NFW halos with embedded spheriodal stellar components on co-planar, prograde orbits in a M31-like 
host galaxy. As the orbit decays due to dynamical friction and the system is disrupted, the stellar 
particles relax to form an extended, exponential disk-like structure that spans the radial range 30- 
200 kpc. The disk scale-length Rd correlates with the initial extent of the stellar component within 
the satellite halo: the more embedded the stars, the smaller the resulting disk scale-length. If the 
progenitors start on circular orbits, the kinematics of the stars that make up the extended disk have 
an average rotational motion that is 30 — 50km/s lower than the host's circular velocity. For dwarf 
galaxies moving on highly eccentric orbits (e ~ 0.7), the stellar debris exhibits a much lower rotational 
velocity. Our results imply that extended galactic disks might be a generic feature of the hierarchical 
formation of spiral galaxies such as M31 and the Milky Way. 

Subject headings: stellar dynamics - methods: N-body simulations- methods: semi-analytical - 
galaxies: kinematics and dynamics - galaxies: halos - galaxies: dwarfs 



1. INTRODUCTION 

Within the context of the hierarchical formation sce- 
nario, galaxies such as M31 are expected to have assim- 
ilated 100-500 smaller mass systems over a Hubble time 
(e.g. Moore et al 1999). Today, stars from these systems 
are expected to be found in the disk, bulge and stellar 
halo (e.g Abadi et al. 2003, 2006, Governato et al. 2004 
and references therein). The fossil signatures of these 
remnants are a veritable treasure trove of information 
about the galaxy assembly process and the nature of the 
merging subunits. 

Detailed investigations of M31 (Ibata et al. 2001, 2004, 

2005, Ferguson et al. 2002, 2005, McConnachie et al. 
2003, 2004, Zucker et al. 2004, Fardal et al. 2006, Font 
et al. 2006, Guhathakurta et al. 2006, Kalirai et al. 

2006, Brown et al. 2006, Chapman et al. 2006) reveal a 
large number of fossil substructures. Intriguingly, many 
of these have disk kinematics and are distributed in a 
gigantic flattened structure surrounding the high-surface 
brightness inner disk (Ibata et al. 2005). The observa- 
tions are consistent with this extended structure being a 
rotating disk-like system whose midplane is coincident 
with that of the inner disk, but which extends > 80 
kpc from the center of M31. The rotational velocities 
of the constituent stars lag the M31 disk rotation by 
~ 30km/s (Ibata et al. 2005). The mass of the extended 
disk is ~ 10% of the inner disk and its stellar surface 
density profile of the structure is consistent with a expo- 
nential decline although there is some uncertainty about 
the value of the radial scale length (cf. Ibata et al 2005; 
Kalirai et al. 2006). Here, we explore whether the ex- 
tended disk could be the remnant of a dwarf galaxy. 

2. CONSTRUCTING OUR GALAXY MODELS 

2.1. The host galaxy 

Our M31 galaxy model consists of three static sub- 
components: a Miyamoto-Nagai (1975) disk, a spherical 
Hernquist (1990) bulge and a spherical Navarro, Frenk 



& White (1997) (NFW) dark matter halo. The corre- 
sponding potentials are: 
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Following Geehan et al. (2005), we set M d = 8.4 x 
10 10 M Q , M b = 3.3 x 1O 1O M , a = 5.4 kpc, b = 0.26 
kpc and c = 0.61 kpc. The M31 halo parameters are 
M & (raoo) = 6.8 x 10 n M Q , r 200 = 180 kpc, r s = 8.18 
kpc, which yields a concentration of C200 = f20o/f s = 22. 

2.2. Generating N-body satellite galaxies in equilibrium 

Our satellite models have two components: (i) a spher- 
ical NFW dark matter halo, (ii) an embedded spherical 
"stellar" component that reproduces a King (1966) den- 
sity distribution and approximates the light profile of the 
dSph galaxies of the Milky Way (MW) and M31 (Irwin 
& Hatzidimitriou 1995, McConnachie & Irwin 2006). 

We use the code of Kazantzidis et al. (2004, 2006) 
to generate NFW satellite halos with 10 6 particles, c = 
r vif/ r s = 20 and an exponential cut-off at the virial ra- 
dius. We calculate the relative energy, e = "J" — v 2 /2, 
and the relative potential, $ = —&nfw + ^nfw(oo), 
for all satellite particles. Next, we generate n bins in en- 
ergy and, for each bin, we select f+(e)/fNFw(c) X N(e) 
"stellar" particles, where N(e) is the number of NFW 
particles within (e, e + de) and de = e max /n. Here, /* 
and Jnfw are the isotropic distribution function for the 
stellar and the dark matter particles calculated from 
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Fig. 1. — Density profile of the stellar component normalized to 
the central censity (po). Solid and dotted lines show respectively 
the profiles at t = and after 14 Gyr of evolution in isolation (i.e 
in absence of external forces). Our method to embed a "stellar" 
profile within a NFW halo produces systems in equilibrium. 



or a King profile. The result is AT* dark matter particles 
with a density distribution 
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r c and rt are the core and tidal radii, K is an arbitrary 
constant, and we have fixed r t /r c = 5 (sec Bullock & 
Johnston (2005) for more details). By construction, the 
stellar component does not contribute to the potential of 
the satellite galaxy, nor does it influence its evolution. 

2.3. Satellite structural and orbital parameters 

We consider satellites with mass M[r v i r ] = 5 x 1O 9 M 
(model HI), the estimated mass of typical Local Group 
dwarf galaxies, and M[r vir ] = 5x 10 10 M Q (model H2), 
at the top end of minor mergers in merger trees for a 
M31-like galaxy. Different King profiles are used for 
the tracer stellar particles, characterized by the quantity 
r c /r s . This determines the compactness of the stellar 
component in the halos. We select r c /r s = 1.0,0.5,0.1 
(Kl, K2 and K3, respectively). The number of stellar 
particles available to trace the King profiles decreases as 
r c /r s decreases, so that A 7 * = 3.1 xlO 5 , 1.3xl0 5 , 1.5xl0 4 
for Kl, K2 and K3, respectively. 

Satellites are placed in the host galaxy at an initial 
apocenter of r a = 75 kpc. The velocity is chosen to 
obtain the desired orbital eccentricity. We explore (i) an 
initially circular orbit, and (ii) a highly eccentric orbit 
(e = [r a - r p ]/[r a + r p ] = 0.71, or r a /r p = 1/6, where 
r a and r p are the apo and peri center distances). The 
latter corresponds to the typical eccentricity of satellite 
galaxies in cosmological simulations (Ghigna et al. 1998). 
For this Letter, we only study satellites whose orbits are 
co-planar and prograde with the M31 disk. 

3. N-BODY SIMULATIONS 



We refer the reader to Peharrubia et al. (2006) for a 
detailed description of our N-body code and the numer- 
ical parameters and techniques used in our simulations. 
Fig-fflshows the density profiles of Kl, K2, K3 at t = 
(solid lines) and after evolving them for 14 Gyr in isola- 
tion. This shows that (i) the method outlined in H2.2l pro- 
duces N-body systems in perfect equilibrium (to within 
Poisson fluctuations) and (ii) the selection of the N-body 
code's numerical parameters is appropriate. 

4. RESULTS 
4.1. Disruption of satellite galaxies 

The process of disruption occurs as the satellite galaxy 
sinks to the central regions of the host galaxy (owing to 
dynamical friction) and tidal forces strip the dark and 
stellar particles. The sequence for model K2H2 can be 
seen in Fig. |5J the left column shows the surface den- 
sity profile of stars (filled circles) and dark matter (open 
circles) at different times. The right column shows the 
projection of stellar (black dots) and dark matter (grey 
dots) particles onto the disk plane of the host galaxy. 
The top panel shows the satellite galaxy at a time when 
most of the stellar material is still bound (see the corre- 
sponding right-hand panel), so that the surface density 
profile is fairly peaked at the satellite location. As time 
passes, the stellar particles start to form tidal streams, 
which can be seen as bumps in £(ii). As the debris con- 
figuration relaxes (mixes in phase-space), the curve £(i?) 
smooths and approaches an exponential-like surface den- 
sity profile over the interval 30 < R < 200 kpc. The solid 
line in the left panels shows the profile of an exponential 
disk with £ = So exp(— R/Rd) and Rd = 30 kpc. Since 
the stellar component is initially tightly bound, then the 
stellar surface density distribution is considerably steeper 
than that of dark matter particles. 

The time-scale for the extended disk to form is a com- 
bination of the decay and relaxation times (idee and t te \ , 
respectively). For massive satellites (eg. Fig. the 
satellite sinks and disrupts within <dec — 3 Gyr. Thus, 
the relaxation and decay times are only comparable in 
the inner regions of the host galaxy. Debris at R > 120 
kpc are not relaxed even after a Hubble time. 

4.2. Spatial distribution of debris 

While the above results are for a specific case, the 
emergence of an extended exponential disk appears to 
be a generic feature of the disruption of satellites on co- 
planar orbits. In Fig.|3 we show the surface density pro- 
files of the stellar and dark matter particles after 14 Gyr 
of evolution for various combinations of satellite masses, 
embedded stellar systems and orbital eccentricities. The 
left and right columns show satellite galaxies on circular 
and highly eccentric orbits, respectively. Top and bot- 
tom panels are for satellite galaxies with initial masses 
of 5 x 10 9 M Q and 5 x 10 10 M Q , respectively. 

In the case of a lower-mass satellite on an initially cir- 
cular orbit, the surface density of stellar debris is not 
smooth even after 14 Gyr and does not resemble expo- 
nential. This is because the satellite has a long orbital 
decay time. The system has disrupted relatively recently 
and the debris has not yet relaxed. This is in contrast to 
the case of a massive satellite galaxy on a highly eccentric 
orbit (left-bottom panel). The dynamical friction is so 
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Fig. 2. — Left column: Surface density as a function of radius 
for the model K2H2 on an initially circular orbit. Stellar and dark 
particles are shown with filled and open circles, respectively. Each 
panel corresponds to a different time. The solid lines shows the 
surface density of an exponential disk S = Eo exp(— R/ Rd) with 
Rd = 30 kpc. Right column: Distribution of stellar (black dots) 
and dark matter (grey dots) debris at the same snap-shots. Note 
how remarkably the stellar debris resemble an exponential disk 
after the satellite galaxy has been totally disrupted. 



Fig. 3. — Surface density of debris as a function of radius after 
14 Gyr of evolution. Models Kl, K2 and K3 are represented by 
filled circles, open circles and stars, respectively. Left and right 
columns show the distribution of debris of a galaxy on a circular 
and a highly eccentric orbit, respectively. Strong solid lines show 
the distribution of dark matter particles for comparison. Note that 
the resulting stellar density depends on both, the dwarf orbit and 
its initial distribution within the halo. 



efficient that the satellite sinks to the galaxy center and 
disrupts within ~ 2 Gyr. Since there is no difference in 
the disruption time-scales for the dark and stellar par- 
ticles, both distributions are practically the same: an 
exponential with a large radial scale Rd — 105 kpc, inde- 
pendent of how compact the stellar system was initially. 

The intermediate cases in the bottom-left (massive 
satellite, circular orbit) and the top-right (average-mass 
dwarf, eccentric orbit) panels show exponential stellar 
profiles over finite radius: 30 — 120 kpc for the circular 
orbit, and 60— 150 kpc for the eccentric one. We find that 
the more compact the stellar component, the steeper the 
final stellar surface density profile. The scale-lengths of 
the exponential profiles are Rd — 50,30 and 7 kpc for 
models with initial r c /r s — 1.0,0.5 and 0.1, respectively, 
independent of initial orbital eccentricity. 

4.3. Debris kinematics 

The extended disk in M31 was primarily identified by 
its rotating kinematic signature. In Fig. 0] we plot the 
average radial and azimuthal velocities of the satellite 
debris for those simulations that form an exponential- 
like disk after 14 Gyr (left column: massive satellite on 
a circular orbit; right column: typical-mass satellite on 
an eccentric orbit). For comparison, we also show the 
circular velocity curve of the host galaxy (dashed lines). 

From the panels showing the mean galactocentric ra- 
dial velocities, we see that the remnant of both satellites 
are practically relaxed, in that the stream motion in the 
radial direction is small. More interestingly, the extended 
debris exhibits a net rotation about the host galaxy cen- 
ter, independent of the satellite's initial orbital eccentric- 
ity. For the satellite galaxy on an initially circular orbit, 
the rotational velocity of the debris is ~ 20 — 50km/s 
lower than the host circular velocity for R > 30 kpc. 
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Fig. 4. — Average azimuthal (top panels) and radial (bottom 
panels) stellar velocities as a function of radius for models Kl, 
K2, K3 (same notation as in Fig. |3J. The left and right columns 
show satellites moving on a circular and a highly eccentric orbit, 
respectively. Note that the velocities of stellar debris are insensitive 
to their initial distribution within the dwarf halo but they reflect 
the orbital properties of the progenitor galaxy. 



The satellite on an initially eccentric orbit shows a much 
lower rotational velocity, typically about lOOkm/s lower 
than v c for R > 40 kpc. Also, the stellar kinematics 
trace the dark matter debris, independent of the initial 
distribution of stars within the dwarf galaxy. 

5. IMPLICATIONS 

We find that stellar debris, resulting from the disrup- 
tion of dwarf galaxies on prograde orbits which are co- 
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planar with the host galaxy disk, relax into an extended 
exponential disk. The disk scale-length is strongly cor- 
related with the initial distribution of stars inside the 
dwarf halo, such that more embedded stellar profiles have 
smaller "disk" scale-length. For the simulations shown 
here, Rd = 50,30 and 7kpc for the initial stellar com- 
pactness parameter r c /r s = 1.0, 0.5 and 0.1, respectively. 
The stellar debris exhibit rotation, even in the case of 
dwarf galaxies on initially highly eccentric orbits. The 
average rotational velocity is around 30 — 50km/s lower 
than the host circular velocity if the satellite initially 
moves on a circular orbit. For highly eccentric orbits, 
the resulting rotational velocity is a factor ~ 2 lower. 

In the specific case of M31, Ibata et al. (2005) find the 
extended disk of M31 has a rotational velocity that is ~ 
30 km/s lower than v c . Based on our results, the observed 
kinematics favour a scenario where the disk progenitor 
was a massive satellite galaxy on a low-eccentricity, low- 
inclination orbit, possibly similar to the progenitor of the 
Monoceros stream in the MW (Penarrubia et al. 2005). 

How common should such extended disks be? ACDM 
simulations find that spiral galaxies experience several 
mergers of massive, 0.01 — 0.1 Mh os t dwarf galaxies (e.g. 
Gao et al. 2004), many on low eccentricity orbits (Ghigna 
et al. 1998). That extended disks appear aligned with 
the inner disk is not surprising taking into account the 
strong decrease in orbital inclination suffered by mas- 
sive satellites in flattened systems (Quinn & Goodman 
1986, Penarrubia, Kroupa, Boily 2002), as well as or- 
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bital circularization (Jiang & Binney 2000). Thus, in a 
ACDM work frame, our results indicate that extended 
disks might be a common result of the hierarchical for- 
mation of spiral galaxies such as the MW and M31. 

Would accretion destroy the inner disk? Our simu- 
lations assume a static, non-responsive host galaxy po- 
tential. Several works (Toth & Ostriker 1992, Walker, 
Mihos & Hernquist 1996, Velazquez & White 1999) have 
shown that, for sufficiently massive satellites, the inner 
disk would be disrupted by the merger event. To anal- 
ize this possibility we have carried out self-consistent N- 
body simulations of the infalling systems considered in 
this Letter. Our preliminar results indicate that the in- 
ner disk is not strongly affected by these accretion events 
since, owing to the nearly exponential mass-loss rate of 
NFW halos, the dwarf galaxies are disrupted before dy- 
namical frictions brings these systems close to the inner 
disk. The results from these simulations, together with 
the analysis of the response of the inner disk, will be 
shown in a subsequent paper. 
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